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In the first part of this talk, I will summarize recent developments in the study of the chiral-odd spin- 
dependent parton distributions hi(x,Q 2 ) and h^{x,Q 2 ) of the nucleon, in particular, (i) Next-to-leading 
order Q 2 evolution of h\(x,Q 2 ) and (ii) Leading order Q 2 evolution of the twist-3 distribution hi,(x,Q 2 ) 
Q\ • and the universal simplification of the Q 2 evolution of all the twist-3 distributions in the large N c limit. The 

' second part of this talk will be devoted to a systematic analysis on the light-cone distribution amplitudes 

' of vector mesons (p, w, <j>, K* etc) relevant for exclusive processes producing these mesons. In particular, 

twist-3 distribution amplitudes are discussed in detail. 
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1 Introduction 



High energy processes can be classified into two categories, inclusive and exclusive processes. Quark- 
gluon substructures of hadrons involved in these processes reveal themselves as a form of parton dis- 
tribution functions in the inclusive processes, and light-cone distribution amplitudes in the exclusive 
' proceeses. Understanding on both quantities constitutes a crucial step for the QCD description of 

the high energy processes. In this talk, I will summarize our recent studies on the quark distribtution 
functions in the nucleon and the light-cone distribution amplitudes for the light vector mesons. 

Spin dependent parton distribution functions for the nucleon measured by the polarized beams 
and targets represent "spin distributions" carried by quarks and gluons inside the nucleon. They are 
functions of Bjorken's x which represent parton's momentum fraction in the nucleon and a scale Q 2 
0^ ' at which they are measured. Untill now, most data on the nucleon's distribution functions have been 

obtained through the lepton-nucleon deep inelastic scattering (DIS). The chiral-odd distributions, 
/ii,l(x, Q 2 ), are the new type of distribution functions which have not been measured so far: Due 
to the chiral-odd nature, they decouple from the inclusive DIS. They can, however, be measured 
by the nucleon-nucleon polarized Drell-Yan process and semi-inclusive DIS which detect particular 
hadrons in the final state. They will hopefully be measured by planned experiments using polarized 
accelerators at BNL, DESY, CERN and SLAC etctl. In particular, RHIC at BNL is expected to 
provide first data on these distributions. 

In the study of these distribution functions, perturbative QCD plays an important role in predict- 
ing their Q 2 -dependence: Given a distribution function, say h\(x, Qq), at one scale Qq, perturbative 
QCD predicts the shape of h\(x,Q 2 ) at an arbitrary scale Q 2 . This Q 2 evolution is necessary not 
only in extracting low energy hadron properties from high energy experimental data but also in 
testing the x-dependence predicted by a non-perturbative QCD technique or a model with the high 
energy data. In the first part of this talk, I will summarize our recent studies on the Q 2 -dependence 
of h ljL (x,Q 2 ). 

Light-cone distribution amplitudes (wave functions) for the vector mesons (p,uj,<f>, and K*) 
appear in various exclusive processes producing these vector mesons in final states, such as B decay, 
B -> tvV, B -> l + £~V, B -> -fV, and electro-production, e + N ^ N' + V. (Study on the wave 
functions for pseudoscalr mesons is less involved, and has been done by many works.) Analysis on 
the wave functions is indispensable to test applicability of perturbative QCD to exclusive proceeses. 
In particular, test of the standard model through the rare B decay requires the knowledge on these 
wave functions. In the second part of this talk, we present a complete classification of the two-particle 
(quark- antiquark) wave functions for the vector mesons based on twist, chirality and spin. This can 
be done in parallel with that for the nucleon's parton distribution functions. In particular, for the 
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twist-3 wave functions, we identify the contribution from the three-particle (quark-gluon-antiquark) 
twist-3 distribution amplitudes, using QCD equation of motion. The renormalization and the model 
building for the twist-3 wave functions can be/should be done starting from these exact relations, 
which is discussed by Tanaka in the workshop. 



2 Distribution Function of the Nucleon in Inclusive Processes 

2.1 Chiral-Odd Distributions hi^(x,Q 2 ) 

Inclusive hard processes can be generally analyzed in the framework of the QCD factorization the- 
oremcl. This theorem generalizes the idea of the Bjorken-Feynman's "parton model" and allows us 
to include QCD correction in a systematic way. Here I restrict myself to the hard processes with 
the nucleon target, such as deep-inelastic lepton-nucleon scattering (DIS, I + p — > V + X), Drell-Yan 
(p + p' — > + X), semi-inclusive DIS (I + p — > I' + h + X). According to the above theorem, the 
cross section (or the nucleon structure function) for these processes can be factorized into a "soft 
part" and a "hard part" : The soft part represents the parton (quark or gluon) distribution in the 
nucleon and the hard part describes the short distance cross section between the parton and the 
external hard probe which is calculable within perturbation theory. For example, a nucleon struc- 
ture function in DIS can be written as the imaginary part of the virtual photon-nucleon forward 
Compton scattering amplitude. (Fig. 1 (b)) According to the above theorem, in the Bjorken limit, 
i.e. Q 2 , v = P-q — > co with x — Q 2 /2v — finite, (Q 2 = —q 2 is the virtuality of the space-like photon, 
P is the nucleon's four momentum), the structure function can be written as 



W(x 7 Q 2 ) 



E 
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,a s ^ 2 )W(y^ 2 



(1) 



where $ a represents a distribution of parton a in the nucleon and H a describes the short distance 
cross section of the parton a with the virtual photon, /i 2 is the factorization scale. In Fig. 1(b), <I> a is 
identified by the dotted line. (Fig.l (a)). Similarly to DIS, the cross section for the nucleon-nucleon 
Drell-Yan process can also be written in a factorized form at s = (Pa + Pb) 2 , Q 2 — > oo with a fixed 
Q 2 / s (Pa,b are the momenta of the two nucleons, Q is the momentum of the virtual photon): 



da f d y» f *Vt>&* (-,-,Q 2 -%,a s (^ 2 )) $ a (y a ,fJ, 2 )& 

a b Jx a J x b \Va Vb M / 



(2) 



where the two parton distributions, <& a and $ b , 
by dotted lines in Fig. 1(c). 



for the beam and the target appear as was shown 





(a) 



(b) 



(c) 



Figure 1: (a) Quark distribution function, (b) Nucleon struction function in DIS. (c) Cross section for the nucleon- 
nucleon Drell-Yan process. 



As is seen from Figs. 1(b), (c), the parton distribution can be regarded as a parton- nucleon 
forward scattering amplitude shown in Fig. 1 (a) which appear in several different hard processes. 
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In particular, the quark distribution in the nucleon moving in the +63 directioiipCan be written as 
the light-cone Fourier transform of the quark correlation function in the nucleonri 



<I>"(.*V)=7>- / ^e™ p - z (PS\r(0)rr(z)UPS), (3) 



where IPS') denotes the nucleon (mass M) state with momentum and spin S^, and ip a is the quark 
field with flavor a. In (^|), we have suppressed for simplicity the gauge link operator which ensures 
the gauge invariance and \ ll indicates the operator is renormalized at the scale [i 2 . A four vector is 
decomposed into two light-cone components a* = -^(a ±a 3 ) and the transverse component a±. In 

(||), z + — 0, z± — 0, and z 2 = 0. T generically represents 7-matrices, T = 7^, 7^75, 1. <I> a (x, fi 2 ) 
measures the distribution of the parton a to carry the momentum k + = xP + in the nucleon, which 
is independent from particular hard proceeses. 

If one puts r = 7 M ,7 M 75, the chirality of ip and tp becomes the same, namely it defines the 
chiral-even distributions. Likewise, putting T — u^u, 1 defines the chiral-odd disributions. For the 
case of the deep-inelastic scattering (Fig. 1 (b)), the quark line emanating from the target nucleon 
comes back to the original nucleon after passing through the hard interactions. Since the perturbative 
interaction in the standard model preserves the chirality except a tiny quark mass effect, the chirality 
of the two quark lines entering the nucleon in Fig. 1(b) is the same. Hence the DIS can probe only 
the chiral-even quark distributions. On the other hand, in the Drell-Yan process (Fig. 1 (c)), there 
is no correlation in chirality between two quark lines entering each nucleon. Therefore the Drell-Yan 
process probes both chiral even and odd distributions. 

The chiral-odd distributions /i"(cc,/x 2 ), h a L (x,[i 2 ) in our interest are defined by putting Y — 
o>^75 in (§):□ 

^e' Xx (PS\r(0)a^i l5 r(*n)\,\PS) 
Ztt 

,2 \ 



= 2[h a l {x^ I ){S^ Pv - S ±vPfi )/M 



+h a L (x 1 fJ 2 )M(p f _ l n l , ~p v n lJ ){S ■ n) 



+h%(x, fi 2 )M(S ±fl n v - S ± „ra M )] (4) 

where we introduced two light-like vectors p, n (p 2 = n 2 = 0) by the relation P^ = p^ + ^-n^, 
P- n = 1, p~ = n+ = 0. If we write P+ = V, p = ^(1,0,0,1), n = ^(1,0,0,-1). Pisa 
parameter which specifies the Lorentz frame of the system: V — > 00 corresponds to the infinite 
momentum frame, and V — ► Mj \[2 the rest frame of the nucleon. is the transverse component 
of S' 1 defined by = (S ■ n)p» + {S-p)rt l + S^. One can show that $ a defined in (g) has a support 
— 1 < x < 1. If one replaces the quark field tp in (||) by its charge conjugation field Ctp T , it defines 
the anti-quark distribution § a . In particular h\ L 3 (x, /i 2 ) in (Q) are related to their anti-quark 
distribution by hi L 3 {—x, /i 2 ) = —ft™ L 3 (x, /1 2 ). 

$ a appears in a physical cross section in the form of the convolution with a short distance cross 
section in a parton level as is shown in ([l]) and (g) . The cross section can be expanded in powers of 
I as 

M M 2 

a{Q 2 ) ~ A(lnQ 2 ) + -_i?(lnQ 2 ) + ^(InQ 2 ) + ■■■, (5) 

where each coefficient A, B, C receives logarithmic Q 2 -dependence due to the QCD radiative cor- 
rection. In order to see how h\ t i,^ can contribute in the expansion (||), it is convenient to move 
into the infinite momentum frame (V ~ Q — > 00). In this limit the coefficient of fti,i,3 in (m be- 
haves, respectively, as O(Q), O(l), 0(1/Q). Therefore if hi contributes to the A term in (pft, Hl 
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spin 


average longitudinal transverse 


twist- 2 
twist-3 


/i .9i hi 
e h L g T 



Table 1: Clasification of the quark distributions based on spin, twist and chirality. Underlined distributions are 
chiral-odd. Others are chiral-even. 



can contribute at most to the B-term, and can contribute at most to the C-term. In general, 
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when a distribution function contributes to hard processes at most in the order of y — 7= j , the 

distribution is called twist-r. Therefore h%, h^, /13 in (||) is, respectively, twist-2, -3 and -4. nnnn 
Twist-2 distribution hx can be measured through the transversely polar^d_^Drell-YanEfQu'El, 
semi-iiiclusive deep inelastic scatterings which detect piontj, polarized baryonsStrEj, correlated two 
pionsEJ. 

From the discussion above, one sees that it is generally difficult to isolate experimentally higher 
twist (t > 3) distributions in hard proceeses, since they are hidden by the leading twist-2 contribu- 
tion (A term in (|^)). However, this is not the case for hL and gT- In particular spin asymmetries, 
they contribute^to the £?-term in the absence of A-term: gT can be measured in the transversely 
polarized DISLL3, and hL appears in the longitudinal versus transverse spin asymmetry in the po- 
larized nucleon-nucleon Drell-Yan processQ. Therefore the (Revolution of gT and hL can be a new 
test of perturbative QCD beyond the twist-2 level. 

Insertion of other 7- matrices in (|J) defines other distributions. In Table 1, we show the clas- 
sification of the quark distributions up to twist-3u There fx, gx,Ti e is defined, respectively, hy 
T = 7^,7^75,1 in (§)■ A similar classification can also be extended to the gluon distributions^. 
The distribution fx contributes to the spin averaged structure functions Fx t 2(x,Q 2 ) familiar in DIS. 
The helicity distribution gx contributes to the Gx(x,Q 2 ) structure function measured in the longi- 
tudinally polarized DIS. By now there has been much accumulation of experimental data on fx and 
gx, and the data on gx triggered lots of theoretical discussion on tha "origin of the nucleon spin" El. 
The first nonzero data on 92 (= .9t — gi) was also reported in Ref. Ej. 

2.2 Next-to-leading order (NLO) Q 2 -evolution of hx(x,Q 2 ) 

As we saw in the previous section, hx is the third and the final twist-2 quark distribution. It has a 
simple parton model interpretation as can be seen by the Fourier expansion of ift in (Q). It measures 
the probability in the transversely polarized nucleon to find a quark polarized parallel to the nucleon 
spin minus the probability to find it oppositely polarized. Here the transverse polarization refers 
to the eigenstate of the transverse Pauli-Lubahski operator 75$j_. If one replaces the transverse 
polarization by the longitudinal one, it becomes the helicity distribution gx. For nqrjxfilalivistic 
quarks, hi(x,(i 2 ) — gi(x,fi 2 ). A model calculation suggests, hx is the same order .as gi .taOEill 

The Q 2 -evolution of hx is described by the usual DGLAP evolution equatioiO. Because of its 
chiral-odd nature it does not mix with gluon distributions. Therefore the Q 2 -dependence of hx is 
described by the same equation both for sirigj^panil nonsinglet distributions. For fx and g\, the 
NLO Q 2 evolutiop_was derived long time agotSEja'Ell and has been frequently used for the analysis 
of experimentscja. The leading order (LO) Q 2 -evolution for hi has been known for some timetl 
In the recepi literature, the_next-to-leading order (NLO) Q 2 -evolution has also been completed by 
two papers^jEa: Vogelsang^ presented the light-cone gauge calculation for the two-loop splitting 
function of hi in the formalism originally used for ffll. WeEj carried outj-Lhe Feynman gauge 
calculation of the two- loop anomalous dimension following the method of Ref. 113 for fx . The results 
of these calculations in the MS scheme agreed completely. In theJnllowing, I briefly discuss the 
characteristic feature of the NLO Q 2 evolution of hi following RefsEj'tEl. 

Analysis of (|4|) gives the connection between the n-th moment of hi and a tower of twist-2 
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operators: 



M n [hi(n 2 )} 



0, 



dx x n hi(x, /i 2 ) 



2M 



(PS ± \Ofa 2 )\PS A 



S±vipa l,a n a i^ 5 (in ■ D) n ijj, 



(G) 



where S± stands for the transverse polarization and 0^(/i 2 ) indicates the operator is renormal- 
izaed at the scale /j, 2 . The contraction with and S± (recall S± ■ n = 0, n 2 = 0) in (^J) projects 
out the relevant twist-2 contribution from the composite operator. ("Twist" for local composite 
operators is defined as dimension minus spin.) By solving the renormalization group equation for 
0„ , one gets the NLO Q 2 dependence of Ad n [hi(fi 2 )] as 



MnMn 2 )} 



a s (n 2 ) 



7<°'/2/3o 



a s {Q 2 ) - a s (n 2 ) ft 




t4 0, V 


4tt fa 




2/3 J 



where a s (Q 2 ) is the NLO QCD running coupling constant given by 



a s (Q 2 ) 
Air 



A,ln(Q 2 /A 2 ) 



/3ilnln(Q 2 /A 



2M 



/3 2 ln(QVA2 



(7) 



(8) 



with the one-loop and two- loop coefficients of the (3- function (3q = 11 — 2/3iV/ and j3\ = 102— 38/3iV/ 

(iV/ is the number of quark flavor) and the QCD scale parameter A. jn and jjp are the one-loop 
and two-loop coefficients of the anomalous dimension j n for O v n S^ v defined as 



In 



4tt 7 ™ Utt/ 



1: 



(1) 



(9) 



If one sets /3i — * and 7„ — > in (Q), the leading order (LO) Q 2 evolution is obtained. 7„ and 
7„ are obtained, respectively, by calculating the one-loop and two-loop corrections to the two-point 
Green function which imbeds O^Sxv To obtain jn \ calculation of 18 two-loop diagrams is required 
fcynman gauge. Since the expression for 7„ is quite complicated, we refer the readers to 



in th 
Refs. 



Ifff 



for them. 



In order to get a rough idea about the NLO Q dependence of hi, we plotted in Fig. 2 7„ 



h(i) 



(jn for h\) in comparison with jh 9(1> (jn 1 for the nonsinglet fi and gi) for Nf — 3,5. One sees 
from Fig. 2 7^ 1 ^ > jl 9 ^ especially at small n. This suggests that the NLO Q 2 evolution of hi is 
quite different from that of fi and gi in the small x region. The relation jn > jn 9 ^ is in parallel 
with and even more conspicuous than the LO anomalous dimensions which read 



jaw urn 



Mo) 



2c '( i+4 S0' 



7„ /9(0) = 2C, 1 



(n+l)(n + 2) 



n+1 \ 



(10) 



To illustrate the generic feature of the Q 2 evolution, we have applied the obtained Q 2 evolution ta 
a reference distribution for gi and hi. As a reference distribution, we take GRSV gi distributions 
and assume hi(x, /i 2 ) = gi(x, fi 2 ) at a low energy input scale U4 2 -.= 0.23 GeV 2 for LO and /i 2 = 0.34 
GeV 2 for NLO evolution) as is suggested by a nucleon modeloii^l. We then evolve them to Q 2 = 20 
GeV 2 and see how much deviation is produced between them. The result is shown in Fig. 3. As is 
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expected from the anomalous dimension, the drastic difference in the Q 2 evolution between hi r 
is observed in the small x region, and this tendency is more significant for the NLO evolution EZtEB'Erl 
(Although gi for u-quark mixes with the gluon distribution, the same tendency in the difference 
from hi is observed for the nonsinglet distribution.) 



300 



200 - 



100 




Figure 2: The NLO anomalous dimension in comparison with 7n ff ^'- This figure is taken from Ref. 




Figure 3: (a) The LO Q 2 evolution of h± (denoted by Su) and gi (denptpd by Am) for the u-quark. (b) The NLO Q 2 
evolution of hi and gi for the M-quark. This figure is taken from Ref. E3 

In Ref. E3, the Regge asymptotics of hi was studied and the small- a; behavior was predicted to 
be h\{x) ~ constant [x — > 0). On the other hand, the rightmost singularity of and are, 
respectively, located at n — —2 and n = — 1 in the complex n plane, which, respectively, corresponds 
to hi(x) ~ x and h\{x) ~ constant as x — > 0. Therefore inclusion of the NLO effect in the DGLAP 
asymptotics gives consistent behavior at x — > as the Regge asymptotics. This is in contrast to the 
(nonsinglet) /1 and g\ distributions, whose LO and NLO DGLAP asymptotics are the same. 

One of the interesting applications of the obtained NLO Q 2 evolution of hi is the preservation of 
the Soffer's inequality^ 2\h c t(x 7 £) 2 )\ < ff(x, Q 2 )+g^(x,Q 2 ). Although the validity of this inequality 
hinges on schemes beyond LOE3, the NLO Q 2 evolution maintains the inequality at Q 2 > Q\ if it is 
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satisfied at some (lowLscale Qq in suitably defined factorization schemes such as MS and Drell-Yan 
factorization schemes So. 

As was discussed in Sec. 2, a physical cross section is a convolution of a parton distribution and 
a short distance cross section. (See (Q) and (||)) For the double transverse spin asymmetry (A^t) 
in the Drell-Yan process, the NLO short distance cross section has been calculated in Ref. O in the 
MS scheme. The analysis on Att combined with the NLO tranversity distribution predicts modest 
but not negligible NLO effectEa 

2.3 Q 2 -evolution of Iil(x,Q 2 ) and its N c — > oo limit 

In general, higher twist (r > 3) distributions represent quark-gluon correlation in the nucleon. Using 
the QCD equation of motion (see ( |33|) later) , one obtains from dJ) the following relation [m q — 0)LJ: 



h L i.r./r) = 2.r j ^|// ;://.//- ) • /' /. ; .r. //" ). (11) 



00 dz"_ ri 
x (PS\\ \iP{uz)ij 5 o- lia gG u a (tz)z»z v iP{-uz)\PS l] ) ) (12) 



I, L [.r.,, 2 ) = ^ f°° %-t ->■''■'■* j vdu j Ids 



z + = and S\\ stands for the longitudinal polarization for the nucleon (S^ — 
S| l = — ^|-n M ). This equation means that consists of the twist-2 contribution and which 
represents quark-gluon correlation in the nucleon. We call the latter contribution "purely twist-3" 
contribution. (Expansion of ( [l2| ) produces twist-3 local operators. See (|l5|) below.) Equation ( pi] ) 
reminds us of the Wandzura-Wilczek relationEZl for gx' 



g T (x,fi ) = / —gi(y,(i )+g T (x,n ). (13) 

Jx V 

For e and gr, one can write down relations similar to ( |l2|). 

The (Revolution of the first and second terms in (|ll|) is described separately. The evolution 
of Iil is quite complicated. A detailed analysis of ( Jl2] ) leads to the following relation for the rt-th 
moment of h$r. 

[{n+1)/2] / ok \ I 



k=2 

1 



R nk = - [^a Xa n X H 5 (in ■ Df- % igG va n v '(in ■ D) n - k ^ - (k -> n - k + 2)] 



(15) 



We note that the number of independent operators {R n k} (k = 2, • • • , [(n + l)/2]) increases with n. 
In the Q 2 -evolution, the mixing among {R n k} occurs and the renormalization is described by the 
anomalous dimension matrix [^ n {g)]ki for {R n k}- If we put the LO anomalous dimension matrix for 
{R n k} as ["ln{g)]ki — {pt s /2'K)[X v \kl corresponding to (||), the solution to the renormalization group 
equation for {R n k} takes the following matrix form: 

[(n+l)/2] 

(P5|||i?„ fc (Q 2 )|P5||>= \l x " /p "] (PS^R^ 2 )^), (16) 

— L J kl 



1=2 



where L = as [®J . X n for hr was derived in Ref. E3. The Q 2 -evolution for ot and e is also described 
by matrix equation similar to (Eq), and the solution was obtained in RefsO'O for gx and in Ref. 
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for e.c3 As is clear from ( |l4| ) and ( [l6| ) .M n [7iL(Q 2 )] and VkL [/i_l(m 2 )] are not connected by a simple 
equation as in the case for the twist-2 distribution (see Although ( |l6| ) gives complete prediction 

for the Q 2 evolution, it is generally difficult to distinguish contribution from many operators in the 
analysis of experiments. 

In order to get a rough idea on the Q 2 -evolution of Hl, we plotted the eigenvalues of X n in Fig. 
4 (right). For comparison, we also showed in the same figure the LO anomalous dimension 7„ /2 
for hi. (Note the differene in convention between (^) and (|l6|).) As is clear from this figure, the Q 2 
evolution of fiL is much faster pthaiU-hat of h\. (See discussion below.) 

It has been shown in RefsaL3eil that at large N c (the number of colors), a great simplification 
occurs in the (Revolution of the twist-3 distributions. Recall X n in Jl6| ) is a function of two Casimir 

operators Cq = N c and Cp 
to 



2N C 



If one takes N c — > oo, i.e. Cp — > N c /2, (14) and (16) is reduced 



M n [h L (Q 2 )} = tfn/0° Mn [h L (fi 2 )], 



In 



2N r 



2(n+l) 



(17) 
(18) 



This evolution equation is just like those for the twist-2 distributions (see (^|)). In Fig. 4 (left), we 
showed the distribution of the eigenvalues of X n obtained numerically at N c — * oo. The solid line 
is the analytic solution in (18), which shows ( |l8| ) corresponds to the lowest eigenvalues at N c — > oo. 
Since ( p7| ) was obtained by a mere replacement Cf — > N c /2 in (|l6|), the correction to the result is 
of 0(l/N 2 ) ~ 10 % level, which gives enough accuaracy for practical applications. 



leading Nc 



exact 




Kgure 4: (Right) Complete spectrum of the eigenvalues of the anomalous dimension matrix for obtained in Ref. 
M . The symbol o denotes the one- loop anomalous dimension for hi. The solid line is the anomalous dimension ( pi) ) 
at large n. (Left) Spectrum of the eige nva lues of the anomalous dimensioiwnatrix for at large N c . The solid line 
denotes the analytic solution given in (ha) . This figure is taken from Ref. P 9 I . 

This large-A^c simplification is a consequence of the fact that the coefficients of R n k in ( |l4| ) 
constitutes the left eigenvector of X n corresponding to the eigenvalue 7^ in this limit: 



[(n+l)/2] 

E 

k=2 



2k 



21 



(19) 
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which implies that all the right eigenvectors of X n except the one corresponding to 7^ are orthogonal 
to the vector consisting of (l — j^y^j- This leads to (|l7|). 

This large- N c simplification of the Q 2 evolution was proved for the nonsinglet 7jt in Ref@ and 
for /i£ and e in Ref.EEI. The corresponding anomalous dimensions for gx and e are, respectively, 



j 4 2(n + 1) 

j=i ■> <■ 1 



f=J J 4 2(n + 1) 



Corresponding to three twist-3 distributions in table 1, there are three independent twist-3 
fragmentation functions.u (Their number is doubled to 6 if one includes final state interactions. 
See Refrl) It has been shown in Ref.c3 that at large N c the Q 2 evolution of all these nonsinglet 
fragmentation functions is also described by a simple evolution equation similar to ([l7|). Therefore 
the simplification of the twist-3 evolution equation is universal to all twist-3 nonsinglet distribution 
and fragmentation functions. 

Toillustrate the actual Q 2 evolution of h^, we have applied ( |l7|) to the bag model calculation 
of /il-Ej (Fig. 5) Fig. 5(a) shows the bag calculation of hja. At the bag scale, purely twist-3 
contribution hi, is comparable to the twist-2 contribution. After the Q 2 evolution to Q 2 = 10 GeV 2 , 
h^ is dominated by the twist-2 contribution (Fig. 5(b)). This can be ascribed to two facts: One 
is the large anomalous dimension (|l^) compared with the LO anomalous dimension of hi (o in the 
right figure of Fig. 4). The other is the presence of a node for hL(x,Q 2 ), which is taken as model 
independent due to the constraint -L 1 dxhi J (x 1 Q 2 ) = 00 which is an analogue of ±he Burkhardt- 
Cottingham sum rule for g%{x, Q 2 )lB. A similar calculation was done for gT in Rcf.c3. Using these 
model calculations, the longitudinal-transverse spin asymmetry, Alt, for the polarized Drell-Yan 
process was estimated in RefEl 
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Figure 5: (a) Bag model prediction for Kl- The dashed line represents the twist-2 contribution to Kl- (b) Bag model 
predreiion for Hl evolved to Q 2 = 10 GeV 2 assuming the bag scale is fj, 2 = 0.25 GeV 2 . These figures are taken from 
Ref.H 

Another simplification of the twist-3 evolution occurs at n — > 00 In this limit, all the twist-3 
distributions obey a simple DGLAP equation (JTt|) with a common anomalous dimension which is 
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slightly shifted from @ and @: 




7„ = 4Ci, V--- +iV c . (21) 



This evolution equation satisfies the complete evolution equation to the 0(ln(n)/n) accurac y0. In 
the right figure of Fig. 4, ( pl| ) is shown by the solid line. One sees that it is close to the lowest 
eigenvalues except for small n. Combined with this n — > oo result, the large- N c evolution equation 
in (17) with (jig) and (EOT) for each distribution is valid to 0((1 /N 2 )\n(n) /n) accuracy. 



3 Light-cone Distribution Amplitudes of Vector Mesons in QCD 

In this section we present a systematic analysis on the light-cone distribution amplitudes (wave 
functions )EJ of the vector mesons (p, u>, (/>, K* etc) following our recent workEl These amplitudes 
are relevant for the preasymptotic correction to various exclusive processes producing vector mesons 
in the final states, such as B meson decay, B — > IvV (semi-leptonic), B — > -fV (radiative), and 
the electroproduction, 7* + N — > N' + V, etc. In particular, we show that the classification and 
analysis of the light-cone distribution amplitudes for vector mesons can be done in parallel with that 
of the distribution functions of the nucleon. (Analysis on the light-cone distribution amplitudes for 
pseudo-scalar mesons is simpler. See e.g. RefsEdH.) For definitness, we discuss the p~ meson wave 
functions. Extention to other vector mesons is straightforward. 



3.1 Definition and Classification 

For the p~ -meson moving in the positive +£33 direction, the light-cone wave functions are defined as 

4>{u^ 2 ) = P+ J |_ e »^(0|S(0)rd(z)|>-(P,A)), (22) 

where \p~ (P\)) stands for the p~ -meson (mass m p ) state with the momentum P and the polarization 
vector e^; P 2 = m 2 , e^ 2 = — 1, P ■ = 0. T denotes generic 7 matrices and z~ is the only 
nonzero componet of the space time cordinate z. The variable u in (j){u) represents a fraction of 
"+" -momentum P + carried by d quark and <f> has a support on < u < 1. Here and below the 
gauge link operator [0,z] = Pexp{ig dtz^A^tz)} which retores gauge invariance is suppressed 
for simplicity. The only difference between the wave function ( p^ ) and the distribution functions 
is that the latter is a forward matrix elements while the former is a vacuum-to-meson transition 
amplitude. In order to classify the wave functions (p^), it is convenient to introduce two light- like 
vectors p and n as was done in section 2.1 in (Q) . They satisfy the relations p-n — 1, P M = p^\ rn 2 n p 

and e| 1 A ' ) = (e( A ) • n)p^ + (eS x ^ ■ p)n^ + We introduce two coupling constants f p and fj by the 
relation 

(0|«(0) 7/1 d(0)|p-(PA)> =f P m p e^\ (23) 

and 

(0\u(0)a^d(0)\p-(PX)) = fj (e«P„ - e ( A >P p ) . (24) 

With these definitions the classification of ( p^ ) can be done based on spin, chirality and twist, as 
was the case for the distibution functions in the nucleon. The only difference is (i) is a vector, 
while Sfi for the nucleon is an axial vector, and (ii) the wave function ( |22] ) should be linear in e^ A \ 
since it is a matrix element between the vacuum and the p meson state. 
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Twist 


2 


3 


4 




0(1) 


0(1/Q) 


0(l/Q 2 ) 


e ii 


h 


HI ' 


93 






ffi , Si 


h 



Table 2: Spin, twist and chiral classification of the p meson distribution amplitudes. Underlined ones are chiral-odd. 

The explicit definitions of the chiral-odd p distribution amplitudes are: 
ge^<OKO)^M|>-(P,A)) 

= if J ( e ±lp» ~ eluip^-LK^ 2 ) + (PuTiu ~ Pvn^){e {x) ■ n)m 2 h { ^ (u, p 2 ) 



and 



dh\'\u,ii 2 ) 

du 



(25) 



(26) 



and the chiral-even distribution amplitudes are defined as 

/ 

= f P m p 



^e^<0|«(0)7 M dMU|p-(P,A)) 



■ n)^{u,p 2 ) +e { ?/Z\u^ 2 ) - ■ n)m 2 p g 3 (u, p 2 ) 



(27) 



and 



— e " (0|u(0)7 M 7 S d(r/ri)| jU |p (P, A)) = - - / p — — — J mpe^ "e^p^ ^ (28) 

All the distribution amplitudes = {0|| , 0j_ , nj? , /in , g± , g± , h%, g^} are normalized as J* du <f>(u) — 
1. The appearance of the derivative form for h^' and g^' in (^6|) and ( ]28| ) is consistent with the 
relation (0|u(0)rd(0)|p-(P, A)) = for T = 1, 7^75. The various factors in front of cj)(u) in (§f)-(§|) 
can be derived by the normalization condition for <p and the QCD equation of motion. Dimensional 
analysis of J25|)-(|2lj|) determines the twist of each wave function, following the same argumnet in 
section 2.1. Table 2 summarizes the classification of 0's. If one ignore the mass difference between 
the u and d quarks, the G-parity invariance leads to <p(u) = <fi(l — u). 



Similarly to the case of the twist-3 distribution functions (see (11) and fll^)), the twist-3 wave 



functions njf , g± and g± contain the higher Fock component describing the multi-particle 
component in the wave function. In order to reveal this, we define three-particle twist-3 quark- 
antiquark-gluon distribution amplitudes as 

2w J 2tt * p (29) 

p. f P e ^ e ^(0\u((n) lal5 gG^(0)d(vn)\p-(P,\)) = Pa\p v e% - P^l]f£ p A{a d , a u ) + . . . 
1 77 (30) 

e^ e *^{0\U(Cn)a afi gG MV (0)d(rin)\p-(P,\)) = 



2n J 2n 
dq f d£ 



2tt / 2tt 



= ^( e(A) ' ri)\PaPn9/3 v -PPPhQL -PaPu9p^ + VpPufft^fIp m P T { a ^ a n) + ■ ■ • » ( 31 ) 
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where g^ v = — p^n^—p^n^, the ellipses stand for Lorentz structures of twist higher than threeE3, 
and V, A, T refers in an obvious way to the vector, axial-vector and tensor distributions, ay, a u and 
\ — oid — a u are the momentum fractions of d quark, u quark and gluon, respectively, and V, A, T have 
a support on a u> a > 0, a u + ad < 1- When the quark masses are the same, the G-parity invariance 
implies that that the function A is symmetric and the functions V and T are antisymmetric under 
the interchange a u «-> ad- This motivates us to define the coupling constants fY p , /jd, f$ p by the 
normalization conditions for {V,A,T}: 

1 rl—oti 

da d I da u (a d - a u ) V(a<j, a u ) = 1, 

.1 /.l-a 



JO 



/ da d / rfa u A(a<j,a„) = 1, 
Jo Jo 



1 fl-a d 

da d / da„ (a d - a u ) T(a d , a u ) = 1. (32) 
o Jo 

Distribution amplitudes defined in this section contain only the light-cone momentum fraction. 
One might ask what about the intrinsic transverse momentum of quarks which causes power correc- 
tion to an exclusive process. Are they missing in this formalism? The answer is no. One can show 
that in calculating the amplitude for an exclusive process the effect of intrinsic transverse momen- 
tum of quarks is systematically and Gauge invariantly incorporated by considering the multi-particle 
(three or more) distribution amplitudes, i.e., it can be reexpressed in the form of higher twist dis- 
tribution amplitudes. 

3.2 QCD Equation of Motion 

The two- and three-particle distribution amplitudes introduced in the previous subsection are not 
independent among one another, but are constrained by the QCD equation of motion. This can be 
most conveniently done by using the identities among the nonlocal operators defining each amplitude: 
For example, relevant identities for the chiral-odd distribution amplitudes are 



{u(x)a f _ iu x l/ d(-x)} = 



l 



d_ 

= i J dvv u(x)x a a a px p 'gG fl 0(vx)d(—x) 

— ix^dp {u(x)d(—x)} — (m u — md)u(x) d(—x), (33) 
u[x)d(-x) - u(Q)d(0) 

I dt dvu{tx)x a a a px >I 'gG fl fj{vx)d(—tx) 
Jo J-t 

+ i f dtd a {u{tx)a aP x p d{-tx)} 
Jo 

+i(m u + m d ) / dt u(tx) $ d(—tx) . (34) 
Jo 

Here we introduced a shorthand notation for the derivative over the total translation: 

(35) 



d 

d a {u(tx)T[tx, — tx]d(— tx)} = -g-^ {u(tx + y)T[tx + y, —tx + y]d(—tx + y)} 



with the generic Dirac matrix structure T. Sandwitching these identities by the vacuum and the 
p meson state, we obtain the relation among two- and three- particle distribution amplitudes. We 
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finally obtain for the chiral-odd distribution amplitudes as 

u 1 

(l-5 + )h^\u) = u Jdv^<P(v)+u fdv~$(v), (36) 

u 

and 

h \\\ u ) = ^ fdv~$(v)- Jdv^$(v) +<y + ^|(u) 

\0 u J 

u u 

+ (zp~T da d da u - T(a d ,a u ), (37) 

au J J 1 - a u - ay 
o o 

where 

U U 

d f f 1 ( d d \ 

+ C3L3- da d da u - I «d ha u - 1 ] T(ay> a«) (38) 

p a« J J 1 - a u - a d \ da d da u J 



with 



5± = JJL'^±^L, CL = ^. (39) 



Here and below we use the shorthand notation u = 1 — u and £ = u— (1 — u) = 2u — 1. According to 
the various "source" terms on the right-hand side of (|36j) and (|37]), one can decompose the solution 
in an obvious way into three pieces as 

hf(u) = hf )WW (u) + hf }9 (u) + hf )m (u), 

h\°\u) = h\* )WW (u) + h\° }9 {u) + h\° )m (u), (40) 

where W| (u) and /ly WW (u) denote the "Wandzura-Wilczek" type contributions of twist 2 

operators, hnP 9 (u) and 3 (u) stand for contributions of three-particle distribution amplitudes and 

hy (u) and h^ m (u) are due to the quark mass corrections. This decomposition can be compared 
with the relation (O) for the twist-3 distribution. In particular, we get 



h\° )WW (u) = 2 I uj dv^^ +uj dv^^j . (41) 

Similarly we obtain for the chiral-even distribution amplitudes as 

u 1 

(l-5 + ^g ( f\u)=u Jdv-^(v)+u Jdv^V(v), (42) 

u 
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and 



9±W = 4 



dv - %>(v) + dv - y>(v 



+ S + (f)j_(u) 



(z p \da d \da u 



1 



^ 3p du l dad l dau 



1 - a d - a u 
V(a d , a u ) 



d 
da d 



1- a d - a v 



- — A(a d , a u ) 
aa„ 



(43) 



where 



*(«) = 20|| («) + <y+f0j.(«) +8-<j>' ± {u) 

1 



da d I da 



da d I da 



1 - aa - a u 
1 

1 — a d — a u 



Q'rf 



Q'rf 



da a 
d 

da d 



a v 



da u 

d 
da u 



V(a d , a u ) 
A(a d , a u ), 



(44) 



with 



f P s _ f P m u ± m d 



f2 



fp 



r V,A 
^3p 



f V,A 

Hp 

f P mp' 



Eqs. ( [42| ) and ( |43| ) again allow the decomposition of g±\u) and gf'(u) into several terms according 
to the source terms: 



(45) 



<#>(„) = - („) + g^(u) + g^ m [u), 

9 { l\u) = g ( l )WW {u) + + g^ m (u), 



(v)WW 



Mm, 



(46) 



where g^ (u) and gj° (u) denotes the contribution from the twist 2 distribution amplitudes 
(Wandzura-Wilczek part), g± 9 (u) and g^ 9 {u) are the contribution from the three-particle distri- 
bution amplitudes V and A. 

The relations (^), ([57|), ( [42] ) and ([l3|) are exact in QCD and form a basis for the rcnormalization 
and model buildings for the twist-3 wave functions. Wandzura-Wilczek parts in (^) and ( [46| ) do 
not mix with other sources under renormalization. Rcnormalization and the model building based 



on the QCD sum rule approach has b 
for the distribution amplitudes in Refs 



,f55 



:ried out in the framework of the conformal expansion 



4 Summary 

In the first part of this talk, I discussed the Q 2 evolution of the chiral-odd spin-dependent par- 
ton distributions hi(x,Q 2 ) and h^x.Q 2 ). The NLO Q 2 evolution for the transversity distribution 
hi(x, Q 2 ) was completed in the MS scheme. This means the Q 2 evolution of all the twist-2 distribu- 
tions has been understood in the NLO level. The resulting Q 2 evolution of hi(x, Q 2 ) turned out to 
cause quite different behavior from the helicity distribution 31(2;, Q 2 ) in the small x region. The LO 
Q 2 evolution for the twist-3 distribution 1il(x,Q 2 ) (and e(x,Q 2 )) was completed. Although their 
Q 2 evolution is quite complicated due to the mixing among increasing number of quark-gluon-quark 
operators, it obeys a simple DGLAP equation similar to the twist-2 distibution in the iV c — > 00 limit, 
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as was the case for the Q 2 evolution of the nonsinglct gT(x,Q 2 ) distribution. The same simplifica- 
tion at N c — + oo was also proved for the twist-3 fragmentation functions. Therefore this large- N c 
simplification was proved to be universal for the twist-3 distribution and fragmentation functions. 

In the second part of this talk, I presented a systematic analysis on the light-cone distribution 
amplitudes for the vector mesons relevant for exclusive processes. Classification can be done in 
parallel with that of the nucleon's parton distribution functions. Constraint relations among two- and 
three- particle distribution amplitudes are derived using QCD equation of motion. These relations 
are exact and must be satisfied for model buildings. Renormalization of the two-particle twist-3 
distribution amplitudes can also be performed based on these relations. 
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